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Abstract 

As an application of Gradient Expansion (long-wavelength) approxima- 

tion, we studied the inhomogeneous universe including the gravitational 

wave(GW). For a plane-symmetric cosmological model, we could implement 

^ the 2nd order expansion and analyze the purely non-linear perturbation of 

O 

GW. The result indicates that the non-linear effect plays an important role in 

^ . . . 

anisotropic stage of the universe. We also confirm that the perturbation from 

the linear source term in gradient expansion is consistent with the result of 
linear perturbation. 

O 

in 

as 

^ I. INTRODUCTION 

I 

OX) Recently, there are greatly successful results in the measurement of the universe. The 

COBE's data [1] is one of the most important result because it indicates the temperature 
fluctuation of cosmic microwave background(CMB), which might be produced in the early 
stage of the universe. We can get the informations of the early universe by analyzing the 
fluctuation of CMB. We also have an alternative way of measurement of the early universe: 
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using the gravitational wave (GW). The most advantageous point to use GW is that direct 
observation of early universe becomes possible because GW can propagate over the last 
scattering surface of the light. Concerning the detection of GW, the LIGO project has been 
developing [2]. Our knowledge of the early universe would become much wider from such a 
measurement in near future. 

Taking into account of developments for these observations, theoretical study in the early 
universe based on general relativity becomes to play much important roles. A purpose of our 
study is to investigate the dynamical feature of inhomogeneous universe. Although it is well 
known from the observation that present our universe is quite homogeneous and isotropic, 
we cannot explain completely this fact. According to "Cosmic No-Hair Conjecture" [3], that 
is attractive idea, the present homogeneity was recovered from inhomogeneous spacetime by 
inflation. However, due to some counter-examples [4], we don't believe that this is true 
in general inhomogeneous case and need to seek the condition for the onset of inflation. 
The point that we should pay much attention is its non-linearity included in the Einstein 
equation. 

There are three ways for analyzing the Einstein equation: searching exact solution, 
solving numerically and using approximation. When we want to investigate the feature of 
inhomogeneities in general case as possible, the approximation is an advantageous method 
unless its validity breaks down. "Gradient Expansion (GE)" method is one of such an 
analysis [5]- [12]. As is shown from its name, this is an approximation which treat the 
spatial derivative term as an expansion parameter. The most characteristic property is that 
one can include the non-linear derivative terms in the expansion parameters. This is the 
reason why one can investigate the non-linearity by using this method. 

GE is valid when spatial derivative terms are smaller than other time derivative terms. 
Provided that L is the scale of spatial inhomogeneity and H is the scale of the time variation 
of the 3-metric 7jj, this means 
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where a is scale-factor and H = a/a. Therefore GE method would become a good approxi- 
mation if the scale of inhomogeneity is larger than the Hubble horizon size. 

Concerning works on GE approximation, the 0th order solutions were investigated by 
Tomita [5]. These solutions were obtained by ignoring all spatial derivative in the Einstein 
equation and regarded as the most general ones in the sense that they include correct 
number of arbitrary spatial functions. Furthermore, he accomplished the formulation of 
the next order expansion keeping to include the right number of functions although it is 
not so applicable one. Note that the same formulation of the expansion scheme is recently 
investigated by Deruelle and Langlois [6]. According to the behavior of 0th order solutions, 
they become the form so-called quasi-isotropic expression, such as 



In spite of ignoring some spatial functions, it is also 0th order solution in GE. Adopting this 
quasi-isotropic solution as the 0th order solution, Salopek and Stewart, Comer et al. studied 
about useful formulation for higher order expansion scheme [7] [8]: Salopek and Stewart built 
Hamilton-Jacobi theory for general relativity and solved the generating functional order by 
order. Comer et al. adopted the synchronous time slice and directly solved the Einstein 
equation expanded by spatial gradient. Both approaches are equivalent mathematically. 
Using these schemes, Croudace et al. [9] investigated the dust collapsing universe to clarify 
the relation with Zel'dovich approximation and found that one of the Szekeres solutions 
was reproduced as the 4th order solution of GE. Nambu and Taruya [10] analyzed the 
feature of inhomogeneous inflationary universe driven by minimally coupled inflaton field 
with arbitrary potential. Furthermore, Tomita and Deruelle studied non-linear behavior of 
perfect-fluid in cosmological case [11]. 

To consider the more generic situation, however, there are some remained issues. Con- 
structing the cosmological model with the gravitational wave is one of the most important 
one. Although Salopek constructed the 0th order solution including the gravitational degree 
of freedom [12], no one investigates the 2nd order expansion. Due to the non-linearity of 
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Einstein equation, the 2nd order correction of GE becomes to play an important role in the 
strong gravitational field. For example, colliding plane wave solutions which are the exact 
solutions of Einstein equation representing the collision of GW, shows us a lot of curious 
features of GW [13]. Hence, we expect that application of GE analysis for GW leads us to 
further recognition of the non-linear features. 

However, one must keep in mind that the present expansion schemes can treat only some 
restricted inhomogeneities. It is due to adopting the quasi-isotropic condition as the 0th 
order solutions. In fact, the quasi-isotropic solutions are not general because some of the 
dynamics of gravitational field are neglected. We can show that a straightforward application 
of the expansion schemes cannot treat the GW completely. 

In this paper, in order to clarify non-linear behavior of the GW in cosmological situations, 
we shall try to proceed to the 2nd order gradient expansion and analyze their evolution. 
Although in some restricted situations, we could complete this analysis without loss of 
feature of GW. 

The paper is organized as follows. In Sec. II we make a short review about the expansion 
scheme of GE analysis and consider how to include GW. We construct a plane-symmetric 
cosmological model with GW and apply the GE analysis in Sec. III. The evolution equations 
up to the 2nd order gradient are derived and evaluated. In Sec. IV, we analyze our model 
using gauge- invariant linear perturbation. Comparing both analysis, we investigate the non- 
linear feature of GW. Summary and discussion are described in final section V. 

II. GRADIENT EXPANSION BY HJ METHOD 

In order to investigate the cosmological model with GW, let us describe a short review 
of the expansion scheme based on Hamilton-Jacobi(HJ) method. Concerning the 2nd order 
expansion including GW, we shall clarify the difficulty in the scheme and consider how to 
avoid it. 

HJ theory is one of the most applicable method for solving the various equations appeared 
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in classical mechanics. For general relativity, HJ method is equivalent to solve Hamiltonian 
constraint and momentum constraint in terms of generating functional [7] [12]. GE analysis 
by HJ method is to expand the generating functional as follows: 

5 = + (2.1) 

where upper script denotes the number of spatial derivative. In the case of the cosmological 
system coupled with n-scalar fields 0a = 1? 2, • • • n), we may take the following ansatz for 
the functional forms [7]: 

5(0) = - - / d'xa'HUc), (2.2) 

K J 

5(2)= j d^xa^[j{4>c)R^K^''{4>c)4>A,^4>r]. (A,5,C = l,2,---n), (2.3) 



where a"^ = y<iet{7jj}. Because of their covariant forms with respect to spatial coordinate 
transformation, the momentum constraint is automatically satisfied order by order. We only 
have to pay attention to the Hamiltonian constraint. Solving Hamiltonian constraint order 
by order, we could get the solution of H^K"^^ . In order to have the time evolution of the 
system, we substitute the generating functionals into the evolution equations. Under the 
synchronous time slice (take the gauge = 1, A^,- = 0), they become 

■ _^6S . 5S 1 , , , 

o4>A a d'-fkt 2 

By solving iteratively, one can obtain the solution in arbitrary order of accuracy. 

However, when we consider the GW, the problem arises; the straightforward application 
of this scheme cannot treat GW. This impossibility comes from the ansatz of forms of the 
generating functionals. We can see that the 0th order functional (2.2) has only one dynamical 
degree of freedoiri of gravity, i.e, only conformal mode(a) and no dynamical mode of GW. 
The 0th order evolution equations (2.4) are obtained by the variation of the functional (2.2), 

2 dH 

7.- = i?7.-, ^^ = -,5^- (2-5) 
Their solutions are expressed in the forms 
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(0a) hij {x) , = 0^ (t + to {x) ) , 



(2.6) 



where hij{x) and to(^) are arbitrary functions of the spatial coordinate. Since we ignore ah 
spatial gradient in the 0th order expansion, the solutions (2.C) arc equivalent to homogeneous 
flat FRW except for its spatial dependence. But due to lack of the GW mode in the 0th 
order solution, even though proceeding to the 2nd order, we cannot recover the dynamical 
feature of GW completely. Therefore we must require a bit device for the application of the 
scheme. 

One of the easiest way to avoid such a problem is to consider a situation in which it is 
possible to include the GW by imposing some symmetry. In other words, we look for the 
system with GW which is equivalent to a cosmological model with conformal mode and the 
scalar fields. If it was possible, we could execute the expansion scheme and study about 
non-linear feature of GW. 

III. GRAVITATIONAL WAVE IN PLANE-SYMMETRIC COSMOLOGY 

In previous section, we found that the present expansion scheme cannot be applicable 
directly to GW. In this section, we describe a cosmological model to which GE is applicable 
and implement the expansion up to the 2nd order. As is shown later, this model has two 
type of solutions concerning number of dynamical mode of GW. Here, we shall analyze only 
one solution(type A solution) for brevity. For the other solution(type B), we mention in 
Appendix C. 



The main assumption to apply the scheme is to impose plane-symmetry, that is, there 
is spatial inhomogeneity only in the z-direction. Consider the cosmological model with the 
anisotropic metric 



A. Plane-symmetric cosmology 



ds"^ — —dt^ + a^^ijdx^dx^ 
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v 



7u = ih IjT] • (3.2) 

1 

Note that the dynamical variables a, ^, // are function of [z. i) due to plane-symmetry. When 
we introduce a cosmological constant A in this system, Lagrangian density takes the form 
as 



2k 



-6(- +- , -2A+ 



1 



a" a' 
-4— + 2 - 
a \ a 



1 f ^ + r}'^ 



77- 



(3.3) 



where (') denotes the spatial derivative with respect to z . One can show that the system 
(3.3) has the dynamical mode of GW. By taking a week field hmit |f | <C 1, I?] — 1| 1, the 
variables (^,r/ — 1) represent ®- ®-mode of GW propagating in the 2;-direction, respectively. 
The transversality of GW is guaranteed by the momentum constraint. 

Now, we shall replace the variables (f, rj) as V2«;(0, a). The Lagrangian density is rewrit- 
ten in the form 



Ml 



(3.4) 



2«; 



-6 ( - j + R{a) - 2A 



; R{a) = 



a" fa'Y^ 
-4- +2 - 



a 



a 



C 



M 



2 p - 2k 
This system is equivalent to a conformally flat cosmology with the spacetime metric 

ds"^ — —dt^ + ^ijdx\lx^\ ^fij — a^Sij^ (3-5) 

and a scalar field $a = (0(t, z) , a(t, z)) . So we can use the 2nd order GE to analyze the 
GW system (3.4). 



B. GE analysis 



We are in the position to analyze feature of GW by gradient expansion. We assume the 
generating functionals for the system (3.3) as 
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K J 

5(2) = y d3xa3 [JR + LC;™?'™ + MC;™??'™ + iV7];™7]'™] , (3.7) 



where the variables J^L^M^N depend on (f,?]). We should keep in mind that (.) and R{a) 
represent the covariant derivative and 3-curvature associated with the 3-metric = a^Sij. 
Hamiltonian and inoirientum constraint become 

The evolution equations are given by 

a K 5S ^ 2k . 2k (3 10) 

Recall that momentum constraint (3.9) is satisfied by the ansatz of the generating functional. 
We only have to concentrate on the Hamiltonian constraint. The explicit forms of the 0th 
order and the 2nd order Hamiltonian constraint are shown in Appendix A. 

Our first task is to obtain the 0th order solution. There are two type of the solution: type 
A and type B. The type A solution we shall analyze specifically here could be constructed 
if we assume that if, J, L, M, N depend only on 77 (For type B, see Appendix C). Then we 
observe M identically vanishes. The 0th order Hamiltonian constraint (A3) yields 

if(r?) = ^(r?^/2 + r?-^/2) ; ifo ^ (1 < < +^). (3.11) 

Substituting (3.6) into (3.10) with a help of (Al), we get the evolution equation for the 0th 
order: 

^ = e = 0, ii^-AjfH^. (3.12) 
From (3.11)(3.12), we obtain the 0th order solution: 



a(0) = a{z) [sinh 3Ho{t + to{z))f^^ , 

(3.13) 

7J 



tanii ^^(^ + ^o(-2^jj 



where tQ{z)^ ^(z), ^(z) are integration constants. Behavior of this solution in , 77) space is 
illustrated in Fig.l. The meanings of the integration constants are as follows, represents 
arbitrariness of the time-slice and a means inhomogeneity of the conformal mode. The 
remained function ^ expresses the inhomogeneity of GW itself. They become the source 
terms of the 2nd order evolution equations. As we are interested in inhomogeneity by GW 
itself, we shall not consider and a hereafter. Note that if we drop the spatial dependence in 
this solution, it belongs to Bianchi type I homogeneous cosmological model. Its asymptotic 
behavior becomes 

Hot »C 1 (Kasner-like expansion) • • • a ~ t^^^, C — f ? rj ^ t ^ . 
Hot ^ 1 (de Sitter-like expansion) • • • a ~ e^°*, C — C? ~ 1- 

Now putting a{z) — 1, to{z) — and noting that r] does not contain a spatial function, 
we derive the evolution equations up to the 2nd order from (3.10) and (A2): 

a ba"^ 

C^-^WiC). (3.14) 

When we solve these equations, we replace with 0th order quantity We should 

pay attention that the evolution equations for the 2nd order include the non-linear term 
(^')^ except for the linear term ({")• The structure of the equations shows that the non- 
linear terms affects the evolution of conformal mode a and one of the GW mode, 77. Hence 
two dynamical mode of GW (f , J]) cannot evolve independently and couples each other. We 
can understand such the effect from the schematic picture. Fig. 2. The left-right arrows on 
the trajectory denotes polarization of the GW in (^, ?])-space. When we set its polarization 
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orthogonal to the trajectory at an initial time, as the universe expands, the arrow becomes 
to incline due to the non-linear coupling of source term and it is no longer orthogonal to the 
trajectory. 

The next task is to solve the evolution equations (3.14). Assuming that the contribution 
of the spatial derivative terms is small, we can linearize the evolution equations: 



(3.15; 



where 5a, 5^, 5rj are perturbed quantities induced by the spatial gradient terms. Substituting 
them into (3.14), we get the linearized equations. The solution can be written easily in the 
integral forms: 



Jo 



Sa — dt 
Jo 



a^J] 



(3.16) 



kL 

6a' 



r\2 



where we choose the integration constant so that the perturbation vanishes at t = 0. In order 
to analyze them, we need L. The second equation of the 2nd order Hamiltonian constraint 
(A4) becomes 

1 



Then L is written in the integral form as 



-2 



(3.17) 



L 



73/2 _ „-V3/2^ 3 ("^ dlj^ fn^/"^ _ Ti"^/^ 



2\ 3 



where we set integration constants so that L vanishes at initial singularity t 
We estimate the asymptotic behavior of L at Hot »C 1 and Hot ^ 1 : 



(3.18) 
0(77 -)■ 00). 



L ~ 



"8(v^+l)re_ffo'^ 



oc-t7f+\ (ifot« 1) 



(3.19) 
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where Loo is a numerical constant defined by 



Lnc, — 



1)2/3 ^0.28; ^(^)^ 



25(1) A s7J-5(s + i)| 

Using the result (3.19), asymptotic behavior of perturbation can be evaluated: 

1. H{)t <C 1 (Kasner regime) : 

(){oc+^t(|"), ()// oc +^^+^(^')^ 5a oc (3.20) 

The result shows that all modes of perturbation grow in time. The growth rate of 
5^ is the same as the result of the usual cosmological perturbation (we will see next 
section). However St]^ da grow as a purely non-linear effect because they evolve through 
the non-linear source term by the 0th order inhomogeneity((f')^). (3.20) indicates that 
non-linear effect such as the mode-coupling will play a significant role in the anisotropic 
universe. 

2. Hot ^ 1 (de Sitter regime) : The perturbations continue to grow by the time ifo^* ~ 1- 
When the Kasner regime ends, they reach certain values Sa^,S^^,Sr}^ . Then the 
universe begins de Sitter expansion (77 — )■ 1) and the perturbations become 



5t] ~ 5t]^ H ^ — -{t] 



T ^ * 



(3.21) 

t. 



In this regime, the perturbations cease to grow and approach constant values. The 
perturbations caused by non-linear source term cannot grow due to rapid cosmic ex- 
pansion. 

To confirm this analysis, we also solved eq.(3.16) numerically. The results are shown 
in Fig 3a. The lines in the figures represent the evolution of amplitude of GW in gauge- 
invariant expressions (For the definition of $, ^, see eq.(4.5)). The solid line is perturbation 
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which comes from the hnear source term, the dashed hne from the non-hnear source term. 
We can observe that all amplitude of GW become frozen after they grow by H^t^ ~ 1. 

We comment on the validity of the approximation before closing this section. The validity 
condition can be obtained from evolution equations (3.10). 



5a 

Roughly estimating, they become 



> 



da 



55( 



6rj 



> 



55(2) 



6rj 



1 

^wavelength of GW^ 



(3.22) 



Thus, we recognize GE analysis is valid when the scale of inhomogeneity exceeds the Hubble 
horizon. Moreover eq.(3.22) is always satisfied if it is satisfied at ifg^* ~ 1 (beginning of 



inflation) . 



IV. COMPARISON WITH LINEAR PERTURBATION THEORY 

We saw in Sec. Ill that two dynamical modes of GW (^, r]) does couple through the non- 
linear source term and their coupling plays an important role at the stage of anisotropic 
expansion. In this section, in order to clarify its feature further, we shall analyze the 
linearized GW and compare with the results of GE analysis. Because we want to extract 
the physically meaningful modes of GW in our discussion, we examine in terms of the 
gauge-invariant expression. 

We start with the definition of the gauge-invariant quantities and derive the gauge- 
invariant wave equations. In anisotropic metric 



ds^ ^ -dt^ + b^dz^ + 



^ ^ dx"^ + 2-dxdy H — dy"^ 
rj rj rj 



(4.1) 



we will define the perturbed quantities as follows: 

a ^ ao(t) (1 + 5a(t,^)) , 6 ^ ao(t) (1 + 56(t, ^)) , ^ ^ ^^{t) + 5^{t,z), rj ^ rjo{t) + 5rj{t, z), 

(4.2) 
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where subscript (o) denotes the homogeneous background solutions. Note that the metric 
(4.1) includes the residual variable b which does not appear in (3.1). It is necessary for 
deriving the physically meaningful equations. Consider the gauge transformation x'^ — )■ 
jf^ — -\- (/i = t, z) which preserves plane-symmetry. Under the synchronous time 
slicing, the solution of e'^ can be written containing arbitrarily spatial functions x(-2), f{z) 

£t = -X(2), £z = x'{z)al / dtriQ^ + f{z). (4.3) 



This arbitrariness produces the gauge-dependent modes in the metric perturbation. The 
perturbed quantities are transformed as 

<^e=-eoX, (4.4) 
^fj = -f]oX, 

dt 



Sb^-Hx-X" f%-f/al 



where H — ao/<^o- It is obvious that every quantities containing Sb become gauge-dependent. 
Using this result, we define the gauge-invariant quantities 

$ = (<^^ + ^<^e)-^f^^'^^^^V«, * = (<^e-^<^^)-^ffo-^eoVa. (4.5) 

We derive the gauge-invariant perturbed equations. In the case of type A background 
solution 



6 = 0, rio 



tanh 



(4.6) 



then the gauge-invariant quantities can be reduced as follows: 

$ = 577 - ^Sa, * = S^. (4.7) 
H 

In order to derive the gauge-invariant equations, we use perturbed Einstein equations rep- 
resented in Appendix B. The evolution equation for $ is obtained from (B2) (B4)(B6) and 
^ obtained from (B3): 
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<i+(3if-2^)$+f^yii+ 



?7o \^f)J 



3 /Ho 



2\H 




$ - $'7a^ = 0, 



(4.8) 



^ + {3H - 2—)* - ^"/al = 0. 



(4.9) 



?7o 



Eq.(4.8) and (4.9) are decoupled so that one can solve them independently. In other words, 
a suitable choice of "polarization" does not make the mode-coupling in linear GW. When 
we substitute (3.16) into (4.7), one can easily see that $, ^ in GE analysis no longer evolve 
independently because of the non-linear source term. Therefore we have confirmed that 
mode-coupling of GW that appeared in the analysis by GE can be recognized as a purely 
non-linear effect. 

Now, let us compare the growing mode obtained from both analysis using the gauge- 
invariant quantities As for GE analysis, we may directly substitute the results (3.16) 
into (4.7). When we evaluate the linear perturbation(LP) from (4.8) (4.9), we take the long- 
wavelength limit(A;/ifo ~ 1) in terms of the Fourier expansion. The results in Kasner regime 
are represented in Table I. One can see that the LP contains the same growing mode that 
is generated from the linear source term in GE analysis. When Hot ^ 1, both modes of the 
LP become the same equation as 



We can observe from this equation that the quantities $, \I/ obey the same equation of motion 
as GW in FRW universe, and all growing modes become frozen in the long-wavelength limit. 
This behavior is the same as (3.21). Accordingly the perturbation from the linear source 
term in GE is consistent with the result of LP in the long- wavelength limit. We may also 
recognize this thing numerically (Fig. 3b). The solid, dashed and other lines represent the 
amplitude \I/ of LP for the various wave number k/HQ. After solving numerically with 
the initial condition ^ = 0, we shift the curves so that they start at the origin. In long- 
wavelength limit {k/Ho ^ 1), the LP coincides with the markers, which is the perrturbation 
from the linear source term in GE. 



$ + 3ifo$ = 0. 



(4.10) 
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type A 






LP 




•^V3(^(^0t)'/') ~ do + rf2t^ + ■ ■ ■ 


GE 


+ (non-linear) 


tH^") (linear) 



TABLE L comparison of the growing mode of GW (type A) 

As for type B solution, similar discussion leads us to the same result as type A: Substi- 
tuting the homogeneous solution of (C3) into (4.5), (B2)~(B4) yield 



$ + 3if$ - - $7a^ = 0, (4.11) 
\i' + 3if^ - ^0 (^j^ * - -^"lal = 0. (4.12) 

There is also no evidence of the mode-coupling in linearized GW for type B. 



V. SUMMARY AND DISCUSSION 

In this paper, we have considered the inhomogeneous cosmological model including GW 
in terms of GE analysis. Imposing plane-symmetry, dynamical modes of GW can be replaced 
into that of scalar fields, hence it has become possible to analyze the feature of GW by using 
the present expansion scheme. The result show us that a non-linearity is incorporated as an 
effect of mode-coupling and it plays an important role at the stage of anisotropic expansion. 
We have also confirmed that the GW generated from the linear source term in GE is same 
growth as the result of linear perturbation theory. 

Consider how initial inhomogeneities affects the present universe. During the inflationary 
regime, the scale of the initial inhomogeneity becomes larger than the Hubble horizon. 
After the end of inflation, since the horizon becomes larger, the scale of the inhomogeneity 
reenters the horizon sooner or later. We can observe such the inhomogeneity as fluctuation 
of CMB due to the fluctuation of the gravitational potential (Sachs- Wolfe effect [14]). In the 
presence of polarized of GW initially, our model indicates that the photon of CMB might 
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be polarized through the non-linear evolution. Taking into account of this non-linearity, we 
can in principle get the information of the early universe from the measurement of CMB. 
Because our model is just preliminary one, we need to study the more realistic model further 
to compare with the observation precisely. 

In fact, we should remark that our treating model belongs to the particular case of Bianchi 
type I so that GW can propagate only to the z-direction (see Sec. Ill A). Thus, in order to 
obtain more general consequence and develop the theoretical prediction to the observations, 
it is necessary to extend our analysis to more general case. But we have not succeeded 
in applying the present expansion scheme to such a general case because we cannot treat 
GW as matter field in such a case. In general anisotropic case, GW has a gauge-dependent 
mode-coupling even in linear level and we cannot have any suitable gauge to eliminate this 
coupling [15]. We must consider a new formulation of the expansion scheme. In order to 
clarify inhomogeneous cosmology, we will proceed to study about this formulation as a next 
step [16]. 
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APPENDIX A: HAMILTON-JACOBI EQUATIONS AND VARIATIONAL 

FORMULAE 

In this appendix, we shall derive the 0th and 2nd order Hamiltonian constraint . 
First, from the ansatz of the generating functional form (3.6) (3.7), we obtain the following 
variational formulae: 
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2 . 
= a^H, 



5rj 



= a^H, 




5a 
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ss 



(2) 



5a 



;m 

m ' 



(A2) 



where the covariant derivative (.) is also induced by the 3-metric = a^S, 



Using these formulae, Hamiltonian constraint (3.8) can be expanded by spatial gradient. 
The resulting equations are 



(A3) 



H^^) = 



HJ-All'{H^J^ + Hr,Jr,) = 2^ 



1 ^-2 



if(L - 4%) - Aii^H,{L, - M^) + A^fH^L^ = -^^77 



if(A^ - 4J,,) - 4r?2if^(iV^ _ M,) + 4r?2if^iV^ = -^^^ 



1 ^-2 



(A4) 



if(M - 8J^^) + 8r?2(if^L^ + if^iV^) = 
if = Tf{2H^L + if^M) 

ifj^ = rf{2Hr^N + H^M) 

At first glance, One can observe that number of equations exceeds that of unknown variables 
and they seem to be inconsistent. However, two equations of (A4) are trivial due to internal 
relation, there is no contradiction (see Ref. [7]). When we solve them concretely for type 
A or B, we can further proceed to reduce them to much simple forms. As for type A, we 
assume H, J. K, L, M, N depend only on r]. Then M vanishes identically and (A4) can be 
reduced to 4 equations. 
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APPENDIX B: PERTURBED EINSTEIN EQUATIONS 



Here, we write down the perturbed Einstein equations under the metric (4.1) 

J2 I ^2 



^ ^ dx^ + 2-dxdy H — dy^ 

7] T] T] 



Using the perturbed quantities defined in eq.(4.2), the perturbed equations of motion are 
obtained: 



1 



5d + 5b + 3H{5a + Sb) + \ ^qS^ + fioSrj 
Sa" 

6d + H{5Sa + 56) = — , 

5i + {3H - 2^)(5e + - f 2^(577 - 25r]] + io{25a + 
?7o ??o V '/o / 



Sii+ {3H -2^)Sij + 



)Sr] + 2—Si + ijo{2Sli + Sb) 



] Sa" 


(Bl) 










(B2) 


2 ' 


(B3) 








(B4) 







?7o ?7o ^70 

where H — d^jaQ^ gq — [sinh HQt]^^^ . In addition, we have two constraint equations of 
perturbation corresponding to the Hamiltonian and momentum constraint: 



1 



H{2Sh + Sb)-^\^oS^ + ijoSi]- 



no 



drj > = Sa" /gq^ 



(B5) 



-Sa 



(B6) 



APPENDIX C: THE TYPE B SOLUTION 

We shall investigate another solution of the system (3.3), so called type B. The type B 
solution up to the 2nd order can be obtained by the following assumption in the functional 
form : 

H^H{x), J^J{x), L^r]-'^L{x), M^r]-^M{x), N^r]-^N{x)- x = ^/r]. (CI) 

In similar manner for type A, we start to construct the 0th order solution. Using (CI) 
and (A3) lead us to the solution H{x) as 
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Thus we can get the 0th order time evolution: 



^(0) = C(^)tanh0(t,^) 
z)sech0(t, z), 



(C2) 



(C3) 



where 



4){t, z) + 00 



V3 



In 



tanh -{t + to{z)) 



(C4) 



The trajectory of this solution is depicted in Fig.l. Note that a^^^t^z) becomes the same 
solution as type A ( see eq.(3.13)). This means that the asymptotic behavior of type B can 
also be classified in two regime: Kasner regime (ifgi 1) and de Sitter regime (ifgi ^ !)■ 

We derive the evolution equations up to the 2nd order. The 0th order solution (C3) has 
the spatial function C{z) which represents the inhomogeneity of GW. When we put a[z) — 
1, t{)[z) =0, the evolution equations are reduced by using (A4) (after some manipulations): 

kP{x) fC'Y 



= H 



6a2 



C 



(C5) 



where. 



2k 



L -N 
N L 



(C6) 



2P f C'\' f C"^ 

P{x) = x'L + N + xM, L = ^--^i-], N = PJ- 



P{x) obeys the following equation: 



1 



4(1 + x')H,P, -HP^ —{x' + 1). 

4:K 



(C7) 



Eq. (C6) shows us that the linear source term L appears in the diagonal part of the right 
hand side matrix, on the other hand, the non-linear term N in off-diagonal part. Because of 
their time dependencies, it is obvious that two dynamical mode of GW no longer decouple 



19 



all over the time. From the result in Sec.IV, a purely non-linear effect is also incorporated 
as the mode-coupling of GW in type B. 

We evaluate (C5)(C6). Dividing them into the perturbed and unperturbed part (see 
eq.(3.15)) and rewrite them in terms of the gauge-invariant expressions (4.5). The results 
are 

ft P fC'Y 



4«; 

$ — Sj] + x5^ — / dt 

V Ju 



J] Jt, 0^(1 + x"^) \C J ^ 

We set all perturbation vanishes at t = ti^Hoti <^ 1). We can solve them numerically. Their 
behavior is the same as type A qualitatively: they grow in Kasner regime and become frozen 
in de Sitter phase. 

Finally, we comment on the differences between type A and B. For 0th order trajectory in 
(^, 77)-space, type B can be represented by a portion of a circle, while type A draws a straight 
line pallarel to 77-axis (Fig.l). As a portion of circle in large raudius limit approximates a 
line, behavior of the type B solution for the 2nd order could also be explained qualitatively 
from type A. 
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Figure Captions 

1. Trajectories of two types of the solutions in 77) -space: The straight hne is the 
trajectory of type A solution. The curved line which is a portion of circle corresponds 
to type B solution. 

2. Schematic picture of the mode-coupling of GW in ({, r^j-space: (a)type A; (b)type B. 
Left-right arrows on each trajectories denotes the polarization of GW in {^,r])-space. 
Due to the non-linear coupling, the initial arrows orthogonal to the trajectory become 
to incline and no longer orthogonal at the end of Kasner regime (see eq.(3.14)(C6)). 

3. Evolution of gravitational wave for type A solution : Horizontal axis denotes time 
normalized by Hq^ and vertical axis means the amplitude of GW $, ^ defined in 
(4.5). 

(a) • • • Results of the GE analysis: The solid line represents ^ which corresponds 
to perturbation from the linear source term and the dashed line is $ from the 
non-linear source term. Solving L from (3.18) and substituting (3.16) into (4.5), 
both lines are obtained numerically. 

(b) • • • Coiriparison between GE and LP: The markers show evolution from the linear 
source term in GE which is the same solid curve in (a). On the other hand, the 
solid, dashed and other lines are the amplitude of linearized GW ^ for the various 
wave number k/HQ. After solving (4.9) numerically under the initial condition 
^ = 0, we shift the curves so that they start at the origin. 
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Result of the GE analysis (Type A solution) 
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Comparison between GE and LP (Type A solution) 
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